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The Phase Diagrams of the Schwinger and Gross-Neveu Models with 
Wilson Fermions 
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A new method to analytically determine the partition function zeroes of weakly coupled theories on finite-size 
lattices is developed. Applied to the lattice Schwinger model, this reveals the possible absence of a phase transition 
at fixed weak coupling. We show how finite-size scaling techniques on small or moderate lattice sizes may mimic 
the presence of a spurious phase transition. Application of our method to the Gross-Neveu model yields a phase 
diagram consistent with that coming from a saddle point analysis. 
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1. Introduction 

There has recently been considerable discus- 
sion on the phase structure of theories with Wil- 
son fermions. A system of free Wilson fermions 
exhibits a second order phase transition at K = 
l/2d, K being the hopping parameter and d the 
dimensionality. Discussions concern the extent 
to which this phase transition persists when a 
bosonic (e.g. gauge) field is included. In the 
Schwinger model, with inverse gauge coupling 
squared (3 , the expectation is that there is a line 
of phase transitions in the (/3, k) plane extend- 
ing to the strong coupling limit j3 = 1. This 
critical line is also expected to recover massless 
physics and therein lies the importance in deter- 
mining its nature and position 

The phase diagram of the Schwinger model 
has been numerically determined by two groups, 
which, while in rough agreement regarding the 
location of the critical line, differ in their con- 
clusions regarding its quantitive critical proper- 
ties. Using Lee- Yang zeroes, finite-size scaling 
and other techniques the results of || support 
the free boson scenario, where the model lies in 
the same universality class as the Ising model, 
with critical exponents v = 1, a = 0. This is 
not in agreement with M where finite-size scal- 
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ing of Lee- Yang zeroes on larger lattices provides 
evidence for v = 2/3, a = 2/3. 

We analytically determine the phase structure 
in the weakly coupled regime. We find that the 
zeroes in the Schwinger model with fixed weak 
gauge coupling display unusual behaviour, do not 
accumulate on the real axis and may not, in fact, 
lead to a phase transition. The behaviour of these 
zeroes is, however, such as to mimic the appear- 
ance of a phase transition when a finite-size anal- 
ysis is restricted to small or moderate lattices. 

Applying our method to the d — 2 Gross- 
Neveu model yields results consistent with Aoki's 
saddle point analysis j| and this consistency adds 
confidence to our new analytic approach. 

2. Lattice Schwinger Model 

We consider a d = 2 lattice of spacing a and 
N sites in each direction. For the fermion fields, 
we impose antiperiodic (periodic) boundary con- 
ditions in the I- (2-) direction. The Fourier 
transformed fermion fields then have momenta 
= (2-7T /Na)p^ where pi is half-integer andjb is 
integer. The action is Sqed = Sg + where 
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This fermion matrix consists of free and interact- 
ing parts, M< w ) = M(°) +M( int ) with M^> given 
by C/p = 1 in (§). Here C/ M (n) = exp (ie aA |4 (n)) 
is the link variable. The partition function 
is Z = /PC/exp({-5 G -4 M/) })(dctAf( w )}, 
the integration over fermions having been 
performed. The gauge action is Sq = 
PJ2p 1 - \{ u p + u p) > where U P is the prod- 
uct of link variables around a plaquette, (3 — 
l/e§a 2 and (O) is a pure gauge expectation. 

For j3 = oo, the partition function is simply 
proportional to detM^ ). Then the free par- 
tition function can be written in terms of its 
eigenvalues Xa\p) = 1/2k — Y^u=i cosp^a + 

j 1 sin 2 p^a which are either 2-fold 
(f>2 — or —N/2) or 4-fold degenerate. In 
this case, and with rj = 1/2k, the zeroes are 

Va (p) 
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The lowest zero with finite real part in k has 
p = (±1/2,0). Pinching of the positive real fi- 
nite hopping parameter axis occurs at k c = l/2d 
and application of finite-size scaling to the imag- 
inary parts of the zero gives the critical exponent 
v = 1 (a = then follows from hyperscaling) . 
Therefore the free fermion model is in the same 
universality class as the Ising model in two di- 
mensions and describes free bosons. 



3. The Weak Coupling Expansion 

The weak coupling expansion is obtained by ex- 
panding the link variables (n) as a power series 
in eo- The fermion determinant can then be ex- 
panded giving the 'traditional' form for the weak 
coupling expansion (see e.g., ||). This expansion 
is analytic in 7/ with poles at r\ = r\i^ Q. 

An alternative formulation of the partition 
function may be obtained by writing the fermion 
matrix as M^ w ^ = r\ + H. The determinant 
detM^ w ^ = det(?7 + H) is a polynomial in r\ for 
finite lattice size. Therefore its pure gauge ex- 
pectation value is also a polynomial in r\ and may 

be written (detM^) = (r) — 7?i). Here % 

represent r\ a (p) and are the zeroes in the complex 
1/2k plane. We may write a 'multiplicative' weak 



coupling expansion as 
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where Aj = r\i — r?|° are the shifts that occur 
in the zeroes when the gauge field is turned on. 
Note that (||) is analytic in rj with poles at r\\ . 
Since the multiplicative expression (|^) must be 
the same as the 'traditional' weak coupling ex- 
pansion, the residues of the poles must be equal. 

Equating these residues gives the the shifts in 
the positions of the zeroes in the two fold degener- 
ate case (the case of four fold degeneracies is more 
complicated but we expect the lowest zeroes to be 
erstwhile two fold degenerate). In this way, the 
shifts in the positions of the lowest zeroes are de- 
termined to 0(e 2 ). Our calculation of pure gauge 
expectation values is done in Feynman gauge. 

4. Results and Conclusions 
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Figure 1. The would-be phase diagram coming 
from the real parts of the lowest zeroes at N = 24. 

A standard numerical technique for determi- 
nation of the phase diagram is to approximate 
the critical point by the real part of the low- 
est zero for some large lattice. Plotted against 
(3 this approximates the phase diagram. In Fig- 
ure 1 we present such a plot for N = 24 (circles) 
to compare with the results of Q (crosses) also 
at N — 24. The phase diagram of Q for the 
Schwinger model with two fermion flavours com- 
ing from a separate PCAC based analysis is also 
included (squares) for comparison. 
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Figure 2 is a finite-size scaling plot for the low- 
est zero at (3 = 10 and N — 8 - 62 coming from 
our weak coupling analysis (circles). The corre- 
sponding data from the numerical analysis of || 
are also included (crosses). The two lines are lin- 
ear fits to the first and second zeroes for lattice 
sizes 16,20,24 which are those analysed in ||. 
These yield slopes —1.5 and —1.4 respectively, 
corresponding to v « 2/3. A similar plot may 
be made to compare with the numerical results 
of @ who use j3 = 5 and N = 2, 4, 8. There, the 
best linear fit gives v compatable with 1. Fitting 
to our small lattice data also gives v w 1. Thus 
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Figure 2. Finite size scaling of the imaginary 
parts of the first two Lee- Yang zeroes at f3 — 10. 

we find that confining the finite-size scaling anal- 
ysis to the range N = 2 - 8 yields v = 1 in agree- 
ment with Q while a corresponding analysis with 
N = 8-24 gives v — 2/3 in agreement with Q . It 
is clear from the figure, however, that the curve 
does not in fact settle to a finite-size scaling line. 
Instead as N increases, the lowest zeroes cross 
the real axis. The first two zeroes therefore fail 
to accumulate and do not contribute to critical 
behaviour [Q . These zeroes are isolated singular- 
ities, having measure zero in the thermodynamic 
limit. Although the possibility of existence of iso- 
lated singularities and non-accumulation of par- 
tition function zeroes has been known for a long 
time , this is to our knowledge the first instance 
where such behaviour has been observed. 

We have also applied our new analytic tech- 
nique to the two dimensional Gross-Neveu model 



|8J. There we find the expected accumulation of 
zeroes and a weakly coupled phase diagram which 
is consistent with that determined by Aoki using 
a saddle point approach 

In conclusion, we have developed a new method 
to analytically determine the Lee- Yang zeroes of 
weakly coupled theories. In the free case, there is 
a phase transition precipitated by the accumu- 
lation of Lee- Yang zeroes on the real hopping 
parameter axis. In the weakly coupled lattice 
Schwinger model at fixed /3, this accumulation 
no longer occurs for the first couple of zeroes. In- 
stead, the movement of these zeroes for small and 
moderately sized lattices mimics phase transition 
like behaviour. As the lattice size becomes large, 
however, these zeroes move across the real axis, 
and do not give rise to a phase transition. In the 
Gross-Neveu case, our phase diagram is consis- 
tent with Aoki's saddle point analysis ||. 
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